Photocount statistics in mesoscopic optics 
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We report the first observation of the impact of mesoscopic fluctuations on the photocount statis- 
tics of coherent light scattered in a random medium. Poisson photocount distribution of the incident 
light widens and gains additional asymmetry upon transmission through a suspension of small di- 
electric spheres. The effect is only appreciable when the average number n of photocounts becomes 
comparable or larger than the effective dimensionless conductance g of the sample. 
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Since the Anderson's discovery that the propagation 
of a quantum particle can be blocked by disorder |lj and 
subsequent realization that this 'Anderson localization' 
can also take place for electromagnetic waves or pho- 
tons (light) 2], the quest for observing it has become a 
breathtaking adventure 0. [J 0- 0- Although the ob- 
servation of microwave localization in quasi-one dimen- 
sional disordered samples |(J now seems to be accepted by 
the scientific community, the localization of visible light 
in strongly scattering, three-dimensional (3D) semicon- 
ductor powders has been questioned 0. The exten- 
sive subsequent work || has shown that the scattering 
strength of available disordered materials in the optical 
frequency range is not sufficient to clearly distinguish the 
impact of Anderson localization on the commonly mea- 
sured quantities (such as the average transmission coef- 
ficient and the coherent backscattering cone) from the 
impact that would be produced by a weak absorption of 
light in diffuse regime. For this reason, well-established 
experimental results in the field of optical localization in 
3D media are limited to 'precursors' of Anderson local- 
ization that can be observed under conditions of diffuse 
scattering: coherent backscattering and weak localization 
, long-range correlations 0, Q , and universal conduc- 
tance fluctuations 

The mesoscopic optical phenomena that we cited above 
can be understood and discussed in the framework of 
classical physics as of 1905, without appealing to quan- 
tum mechanics. The impact of quantum-mechanical ef- 
fects on the coherent backscattering of light has been 
demonstrated in beautiful experiments on light scatter- 
ing in cold atomic clouds [l0(. The quantum nature of 
scatterers (atoms) had to be taken into account to un- 
derstand the low value of the coherent backscattering 
enhancement factor. Experimental studies revealing the 
quantum nature of light in multiple scattering have been 
reported only very recently ^]J, besides a considerable 
theoretical interest in this subject |l2|. In particular, 
Kindcrmann et al. have predicted that disorder can sub- 
stantially alter photon statistics of degenerate incoherent 
radiation. In the present Letter we report the first exper- 



imental observation of the impact of one of the precursors 
of Anderson localization — mesoscopic, long-range cor- 
relations — on the photon statistics of degenerate coher- 
ent light emitted by a conventional continuous laser. We 
interpret our results in the framework of the semiclassi- 
cal theory of photoelectric detection ^3], which appears 
to be sufficient under conditions of experiments reported 
here. 

As first noted by Einstein 100 years ago the quan- 
tum nature of light can be directly probed by the photo- 
electric effect. As the energy of the electromagnetic wave 
is quantized in portions Hlo (with fi the Planck constant 
and u) the frequency of light), only an integer number n 
of such quanta (photons) can be absorbed by a photo- 
electric effect-based detector during a given time interval 
t. Today's electronic equipment allows us to measure n 
in a wide dynamic range and to determine the probabil- 
ity distribution of photocounts P(n,n), where n is the 
average number of photocounts in the interval t — n/ f 
and / is the average photocount rate (number of counts 
per unit time). P{n,n) carries fundamental information 
about interaction of light with the medium. In a 'random 
laser', for example, P{n,n) can be used to characterize 
different regimes of lasing (lo|. In a different domain 
of physics — mesoscopic electronics — the statistics of 
quasi-particle (electron) counts ('full counting statistics') 
in disordered conductors is also under active study [lij ]. 

In our experiment we measure the distribution of pho- 
tocounts of laser light transmitted through an optically 
dense slab (see Fig. The sample cell (thickness 

L w 0.5 mm) is filled with a charge stabilized aque- 
ous colloidal dispersion of a commercial titanium diox- 
ide powder (Warner Jenkinson Europe Ltd.), particle di- 
ameter ss 200-300 nm, at an initial density of 18 ± 1% 
per volume. To further increase the density we let the 
sealed suspension settle under gravity. Due to the electro- 
static repulsion between the particles the sedimentation 
is asymptotically slowed down and an equilibrium layer 
of approx. 0.2-0.25 mm thickness (volume fraction ca. 
35-40%) is formed after about 10 hours. From diffusing- 
wave spectroscopy measurements we have checked that 
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FIG. 1: Experimental setup: The incident laser beam is fo- 
cused on the vertically oriented sample holder. The glass win- 
dow directed towards the laser is highly absorbing while the 
opposite window is transparent. The windows hold a sealed 
liquid layer (thickness L) of colloidal titanium dioxide sus- 
pended in water. Light transmitted through the sample is 
collected by a light guide and recorded by a single photon 
detector (not shown). 



the particles in this layer remain mobile and undergo 
Brownian motion. Following the approach described in 
reference |9j we estimate the transport mean free path 
at this density to I* = 0.7 ± 0.1 fim. Our sample is 
therefore deep in the multiple scattering regime: a typi- 
cal transmitted photon experiences (L/l*) 2 ~ 10 5 elastic 
scattering events, whereas the coherent incident beam is 
destroyed after a distance I* <§; L and hence does not 
contribute to the measured signal. The glass window di- 
rected towards the laser is highly absorbing (transmission 
coefficient ~ 0.001) in order to suppress multiple reflec- 
tions ]<J while the opposite one is transparent. A fre- 
quency doubled Nd:YV04 laser ('Verdi' from Coherent) 
operating at the wavelength A = 532 nm illuminates the 
sample through a microscope objective that focuses the 
laser beam to a small spot (spot size w 3 fmx) on the 
sample surface. The light transmitted through the cell is 
collected by a light guide (core diameter 5 mm) , recorded 
by a single photon detector with short dead time, and 
processed by a digital photon counter (correlator.com, 
New Jersey, USA). The high temporal resolution (12.5 
ns) of the latter assures that no more than one photon 
is arriving every time step for a typical photon count 
rate / of 9 MHz. In a typical experiment the photon 
trace is recorded over one hour. From the recorded data 
we compute the probability distribution P(n,n) and the 
time averaged correlation function of total transmission 
GAt) = (T(0)T(t))/T 2 -1 following standard procedures 
|l8| . To suppress contributions from slow drifts the data 
is analyzed in thirty second intervals and subsequently 
averaged. Some representative results for different beam 
spot sizes are shown in Figs. Eland EH 

Because the detection process is of probabilistic na- 
ture, the detection of a photon is a random event and 
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FIG. 2: Probability distributions of photon counts for a wide 
incident beam (black squares) and a focused beam (beam 
waist w = 3.4 /xm, red circles) for four different sampling 
times t. The former distribution follows the Poisson law 
(black line), while the latter one is well described by the 
Fourier transform of our Eq. © (red line). Small deviations 
for the shortest r can be explained by the finite detector dead 
time [Tt| . The number of data points for the two longest r 
has been reduced to improve readability. 



P(n,n) is expected to be the Poisson probability distri- 
bution 0|. However, as follows from Fig. [2 this is only 
true when the incident laser beam is sufficiently wide. For 
a focused beam we observe that the distribution widens 
and becomes more asymmetric than one would expect for 
the Poisson distribution. This indicates that additional 
fluctuations exist in the latter case. These additional 
fluctuations are due to the random motion of scatter- 
ers. Since we collect all the transmitted intensity (total 
transmission measurements), and since the surface of our 
sample is much larger than the typical size ~ A of speckle 
spots, one could naively expect spatial speckle to average 
out, as it is indeed the case for large w. However, if w 
is small and the scattering is sufficiently strong, coherent 
interferences of scattered light give rise to weak but long- 
range correlation of distant speckle spots, which acquire 
a synchronous component in their fluctuations. This re- 
sults in enhanced fluctuations of the total transmitted 
signal 0. 

The enhanced fluctuations of the total transmission T 
can be studied assuming that l ight is a classical wave de- 
scribed by Maxwell equations 0, • As long as local- 
ization effects are weak, the statistical distribution Pt (T) 
of T appears to be very close to Gaussian (because the 
many speckle spots contributing to T are only weakly 
correlated and the central limit theorem applies), but 
with enhanced second and non-zero third central mo- 
ments A4 2) and A4 3) , where = ((T — T) k ) /T k . 
Therefore, the characteristic function of T can be ap- 
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FIG. 3: Correlation function of total transmission T for dif- 
ferent beam spot sizes w (distance between 1/e intensity val- 
ues of a focused Gaussian beam). Lines are fits to the data. 



proximated by 

X r(«) = exp (iTq - - ^4V) (1) 

When the measurement of T is not instantaneous but 
involves time integration, the second moment is given by 

El 



r 

M t } = 1 J (l-;)c 2 (t)<« 
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(3) 

The third moment My of the distribution of T can be 
shown to be proportional to the square of the second 
one: M^ 3 ' = aM^ 2 , where the proportionality con- 
stant a = 16/5 for a wide [w S> L) Gaussian beam 
[2fj|. The limit of wide beam has also been studied in 
the previous correlation experiments and the corre- 
sponding correlation function Cg(t) has been analyzed 
theoretically |2l| . In our experiments, on the contrary, 
the beam width w is much smaller than the thickness L 
of the sample (typically, w/L ~ 10~ 2 ). In this situation, 
by performing calculations similar to that of Ref. [2l| we 
find C 2 (t) = (2/3 3 )exp(3t/4i )[l - <f>(y/3t/4t )], where 
$ is the error function. Leaving the discussion of the mi- 
croscopic expressions for g and to for a future publication 
[23, we j us t note here that g scales roughly as 1/w and 
hence the magnitude of the total transmission fluctua- 
tions (~ 1/fl 1 ) can be varied by adjusting the beam spot 
size w. By analogy with the case of disordered waveg- 
uide [||, we will further term g the 'effective' dimension- 
less conductance. The above expression for Czit) with 
g ~ 10 3 and t ~ 10 _5 -10 -4 s provides a good fit to our 
measurements (see Fig. OJ). 

According to the famous Mandel's formula [l3| . 
the statistical distribution of photocounts P(n,n) can 



be obtained by averaging the Poisson distribution 
= TjT) over the distribution Pt(T) of the 
total transmission T, with rj the quantum efficiency of 
the photodetector. The Fourier transform of the Man- 
del's formula with respect to n yields a relation between 
the characteristic functions Xn(o) and xt (<z) • 

Xn(q) = Xt [in(l — e iq )/T\ 



exp 



^-><V- i^V] (3) 
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where the second line is obtained by expanding the first 
one in power series in q, which is justified for large n. 
The second and the third central moments of n in Eq. 
© are 

1 



M (2) = A + M i 2) 



Mi 3) = -k+M. 



(3) 
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The Fano factor F 
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((n 2 ) - n 2 )/n = 1 + nM, 
1, which indicates photon bunching. The photocount 
distribution P(n,n) obtained by the Fourier transform 
of Eq. J2Jl describes our measurements very well (see Fig. 

, when we use the fits to the correlation function C2 (t) 
obtained independently (Fig. to determine M^- 1 and 
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The first terms on the right-hand sides of Eqs. JIJ and 
(0 correspond to the results expected for Poisson dis- 
tribution of photocounts. As we show in Fig. ^ these 
results are recovered for small n -C g, when the photo- 
count distribution is dominated by the shot noise due 
to the discreteness of n and is not sensitive to the ran- 
domness of the scattering medium. Deviations from the 
Poisson-like behavior start to become important when n 
becomes comparable to the effective dimensionless con- 
ductance g. The long-range character of the correlation 
function C2(t) ~ 1/Vr is responsible for new scaling laws 
Mn ~ 1/VB and M„ ~ 1/n in the limit of large n. 
This behavior is well confirmed by our measurements: 
as can be seen in Fig. the data points indeed follow 
the and 1/n asymptotes shown by dashed lines 

[23j |. Mesoscopic fluctuations of the total transmission 
T due to the random motion of scatterers in the disor- 
dered sample become dominant in this regime, whereas 
the shot noise is negligible. We see therefore that the 
transition between small- and large-n regimes in P(n,n) 
is governed by localization effects, the strength of the 
latter being measured by the dimensionless conductance 
9- 

In conclusion, mesoscopic fluctuations of coherent light 
transmission through a random medium produce measur- 
able deviations of photocount distribution P(n,n) from 
Poisson law, provided that the average number of pho- 
tocounts n is comparable or larger than the effective di- 
mensionless conductance g of the random sample. This 
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FIG. 4: Second and third central moments of the photocount 
distribution P(n,n) for the same beam spot sizes as in Fig. El 
Lines are theoretical results £Q and (JKJ. The third moment 
plot is a fit to the data with a = 3.02, 2.59 and 2.00 (for 
curves from top to bottom). 



provides a new tool for studying mesoscopic phenomena 
in random media, a tool that should be particularly valu- 
able in the search for Anderson localization of light. An 
interesting continuation of this work would be to analyze 
P(n,n) for incident light in thermal or in non-classical 
(Fock, squeezed, etc.) states. 
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